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Abstract 

A simple analytical/numerical model has been developed for computing the evolu- 
tion, over periods of up to a few hours, of the current and temperature profile in the 
upper layer of the ocean. The model is based upon conservation laws for heat and 
momentum, and employs an eddy diffusion parameterisation which is dependent on 
both the wind speed and the wind stress applied at the sea surface. Other parame- 
ters such as the bulk— skin surface temperature difference and CO2 flux are deter- 
mined by application of the Molecular Oceanic Boundary Layer Model (MOBLAM) 
of Schlussel and Soloviev. A similar model, for the current profile only, predicts a 
temporary increase in wave breaking intensity and decrease in wave height under 
conditions where the wind speed increases suddenly, such as, for example, during 
gusts and squalls. The model results are compared with measurements from the 
lagrangian Skin Depth Experimental Profiler (SkinDeEP) surface profiling instru- 
ment made during the 1999 MOCE-5 field experiment in the waters around Baja 
California. SkinDeEP made repeated profiles of temperature within the upper few 
metres of the water column. Given that no tuning was performed in the model, 
and that the model does not take account of stratification, the results of the model 
runs are in rather good agreement with the observations. The model may be suit- 
able as an interface between time-independent models of processes very near the 
surface, and larger-scale three-dimensional time-dependent ocean circulation mod- 
els. A straightforward extension of the model should also be suitable for making 
time-dependent computations of gas concentration in the near-surface layer of the 
ocean. 
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1 Introduction 



The temperature of the sea surface, and its relation to that of the adjacent at- 
mosphere and wat er column, are vital ly important para meters for the air-sea 
exchange of heat (IFairall et al. 1996a) and gas species ( McNeil and Merlivat 



1996: lFairall et alJ200(ilWajd et all2004al) . For global estimation of atmosph 



lere- 



ocean heat and gas flux, it is necessary to relate satellite observations of radia- 
tive surface skin temperature T s to the bulk temperature T R of the upper cen - 
timetre or so of the water column ( Schlussel et al. 19901 Donlon et al. 20021 ). 
This bulk— skin temperature difference ATg-s — Tb — T$ is a function of the 
ambient radiative, mechanical, and thermal forcing, and is generally positive 
and equal to a few tenths of a degree, primarily as a result of the cooling of 
the water by outgoing long-wave radiation emitted in a thin layer of the order 
of 1/zm deep. The transition zone from this skin layer to the bulk is controlled 
primarily by molecular heat conduction and has a thickness of the order of 
1 mm. 



A number of diff erent theoretical and numerical models exist for the calcu- 
lation of ATg-s f|Eiflerlll992l Il993t I.Tessup et al]ll99.4 ISoloviev and Schliissell 
1994). These models can be quite complex, since they compute many interre- 
lated parameters, and, since they consider processes with a small time scale, 
are generally in time-independent form. Alternatively, one can use such tech- 
niques as neural network methods t o derive empirically the relation between 
the forcing parameters and /ST B _ S ()Ward and RedfernlE^l . 



In the upper 10-50 m layer of the water column, time- dependent three-dimen- 
sional ocean circulation models may employ simple algebraic parameterisa- 
tions for turbulent fluxes of momentum, heat, and mass. They may also use 



some form of turbulence closure technique (e.g. iBurchard et al. 19981 ). albeit 
at greater computational expense. Even more expensi ye computatio nally are 
direct numerical simulations of vortex structures (e.g. iNagaosa Il999l ). 



In this paper we employ a time-dependent model which should be applica- 
ble to the zone intermediate between the surface millimetre layer where the 
skin effect occurs and the «10m level which can be resolved by oceanic cir- 
culation models. The model is thus useful in describing or parameterising the 
coupling between the different model scales. It uses a simple parameterisation 
of turbulent flux processes, so that it should be numerically stable and com- 
p utationally inexpensive. T he model codes and documentation are specified 



in 



Jenkins and Wardl (|2003f ) 



We compare results of the model with near-surface temperature profiles ob 
tained using the Skin Depth Experimental Profiler (SkinDeEP, IWard et al. 



2004bl ). It is an autonomous profiler, carrying high-resolution temperature 
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sensors to provide a record of the bulk temperature. It was deployed dur- 
ing the Marine Optical Characterisation Experiment (MO CE-5), cruise in the 
Gulf of California during the period 1999 October 1-21 ((Ward and Minnett 
20041 ) . 



2 Model formulation and specification 



2. 1 Physical basis of the model 



The flux of momentum, heat, and mass between the atmosphere and the 
ocean is to a large extent controlled by turbule nce in the atmospheric an d 
oceanic boundary laye rs, although s tratification ( Monin and Obukhoylll954l). 
the Earth's rotation fekman 1905f ). surface waves ((janssenl 19891 : Jenkins! 
1987aL Il989t IWeberlll98.C the presence of laminar l ayers close to the sur face, 



and the short-wave and lone-wave radiation balance (Schlussel etal . 1990) can 
also be important . The effect of turbulence may be modelled by eddy viscos- 
ity (|Madsenlll977l ) or turbulence closure mo dels which have reach ed a great 
degree of complexity and sophistication (e.g. Burchard et al. 19981 ). However, 
one difficulty that may arise when applying turbulence closure models in cal- 
culating properties very close to the sea surface is the choice of boundary 
condition: in general, it may be necessary to specify a roughness length in 
a more-or-less arbitrary fashion, and the values of mean velocity, tempera- 
ture, and concentration near the water surface may have a behaviour which is 
close to singular and thus difficult to resolve numerically. The partial differen- 
tial equations which are obtained from a turbulence closure formulation may 
be rather complex and their solutions may pose analytical and/or numerical 
difficulties. 



The model approach employed in this paper is designed to avoid some of 
the difficulties described above, by combining physical conservation laws with 
some general properties of turbulent boundary layers and empirical obser- 
vations. The physical conservation laws are those of momentum, heat, and 
mass: the rate of change of the momentum of the water column is given by 
the applied wind stress; likewise the atmosphere-ocean heat flux is given by 
(or determines) the rate of change of the total heat content of the water col- 
umn, and the flux of a substance through the sea surface determines the rate 
of change of its concentration in the water column. The empirical observations 
are as follows: the mean velocity (current) at the sea surface is determined, 
not by a relationship between the applied wind stress and the vertical velocity 
gradient, mediated by a viscosity or turbulent eddy viscosity, but, instead, as 
a fraction A ~ 0.02 of the wind speed U. Thus, the value of the wind-induced 
surface current has a relatively stable behaviour, which is in agreement with 
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the practical engineering calculations made in connection with observations 
of the drift of floating objects and of oil slicks, and the measureme nt of near- 
surface currents us i ng moored current m eters and surface drifters (jAudunson 
19791 : lHuanelll979L Il983t I Jenkin j Il987l bh . Although the phys cal reason for 
this behaviour of the surface current may not be strictly evident, we may 
consider that an increasing wind stress will cause an increasing vertical ve- 
locity shear, but it will also cause an increase in turbulent motion which will 
transport momentum more rapidly downwards, thus decreasing the shear. The 
corresponding conditions for temperature and concentration would be equal 
temperatures and concentrations just above and below the water surface, al- 
though corrections need to be applied for the presence of a lam inar boundary 
layer and the thermal radiation balance ( Schfiissel et al. 1990j ). 



The effect of turbulence in transporting momentum, heat, and mass vertically 
can be parameterised in terms of a mixing length proportional to the dis- 
tance from the surface, as in the usual turbulent boundary-layer of Prandtl 
and von Karman. However, we show that if we assume a simple functional 
form (Eq. 1 below) for the time-dependent behaviour of the current, temper- 
ature, etc., combined with the conservation laws for momentum, heat, and 
mass, and the empirical surface boundary conditions discussed above, an ef- 
fective 'mixing length' scale appears as a consequence. Although the physical 
boundary conditions at the sea surface and at the sea bottom are fundamen- 
tally diff erent, there is evidence from observations of wind- in duced surface 
currents (lHuanglll979Lll98 3r I Jenkinsl 119 3 l,Tenkins et aD lf986l. temperature 
(jSoetie and Hub erl ll98nl: iBurchardl l2002l. secti on 7.3.2), and the behaviour of 
air bubbles in the water column ( Thorpe! 1984 ). that a substantial vertical ve- 
locity shear can exist near the surface, consistent with an increase in eff ective 
eddy viscosity with increasing distance from the surface, as in Madsen's 11977 
approach. 



2.2 Predicted current profile 



For an ocean initially at rest, we consider at first the following evolution of 
the current u(z,t) at time t and depth —z: 



u(z,t) = \U exp \Uz/(u 2 t) 



where is the friction velocity within the water column, r = p w u 2 = p a U 2 
being the wind stress. U*, p a , and p w are the friction velocity in the atmosphere, 
the air density and the water density respectively. Equation 1 satisfies the 
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conservation of momentum: 
o 



d 
'dt 



u(z, t) dz = p a U 2 , 



(2) 



where for simplicity we neglect the Coriolis force. 
The current u(z,t) in (1) obeys the equation 

u 2 z d 2 u 



du 
dt 



XU dz 2 



(3) 



with —u 2 zj (XU) playing a role similar to that of an eddy viscosity, as in 
turbulent boundary layer flow near a rigid wall. We may also consider this 
'quasi eddy viscosity' to be a product of a turbulent velocity scale equal to 
and an effective turbulent roughness length equal to —u*z/(XU). 

If the wind speed is allowed to vary with time we can alter (1) so that it 
becomes 



u(z,t)=J X—exp 



XU(t') z 



dt'. 



(4) 



In this case, the partial differential equation (3) does not hold, as the system 
retains a 'memory' of the wind forcing at previous points in time. The memory 
is effectively longer at large depths than near the surface. 



Figure 1, from Jenkins! ( 2001 ). shows the model predictions for the current 
induced by a wind of 20 m/s blowing for 2 hours followed by an increased 
wind of 30 m/s. The current profiles are calculated every 5 minutes. We assume 
that A = 0.02, and = {C D p a /p w ) 1/2 U, with p jp„, = 1/800, C D = A + BU, 
A = 0.8 x 10" 3 , and B = 0.065 x lO^m^s (|Willl982h . 



Effect of the modelled current on surface waves 



Wave spectra have been calculated by using the lDonelan et al. (1985) formu- 
lation for limited fetch, transformed by means of the wave group velocity to 
apply to winds blowing for limited times. Calculations have also been made of 
the effect of the computed vertical current s hear on the wave height req uired 
for wave breaking, according to the theory of Banner and Phillips (1974). Re- 



sults of these calculations (Figure 2) indicate that a rapid increase in wind 
speed will tend to increase the intensity of wave breaking so that the wave 
height will decrease by a small factor, of order 6% according to linear wave 
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u / m s A -1 



1. Exam ple of the evolution of a current profile using the present model. (From 
nkins 2001). © 2001. ASCE. Reproduced by permission of the publisher. ASCE. 
http : //www . pubs . asce . org 

theory, but this factor may be as much as 25% accordi ng to time-dependent 
nonlinear numerical simulations by Banner et ah ( 2000[ ). 



2. 3 Heat flux 



If we now consider the diffusion of heat within the water column, we obtain a 
corresponding formula for the evolution of the temperature T: 



TV A T a. f dT (°~^) 



XU{t') z 



K(f)r (*-*') 



dt', 



(5) 



where T is the initial temperature, and T(0~,t) is the bulk temperature just 
below the surface skin layer. Here we assume that turbulence diffuses heat 
at the same rate as it diffuses momentum. We may adjust the relative values 
of the turbulent diffusion coefficients for heat and momentum by altering the 
value of the parameter A. The model is assumed to be horizontally homoge- 
neous, and the effect of density stratification on the turbulent diffusion of heat 
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Hz 



Fig. 2. Wave spectra from the Donelan/Hamilton/Hui model, an d com- 
puted reduction in amplitude of maximum non-breaking waves. (From Jenkins! 
2001). © 2001, ASCE. Reproduced by permission of the publisher, ASCE. 
http : //www . pubs . asce . org 



is neglected. 
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2.4 Characteristic behaviour of the model 



It may be considered that the model presented in this paper produces results 
which are at variance with the turbulence closure and eddy diffusion models 
normally employed. In particular, the evolution of the current, temperature, 
or concentration of a tracer, driven by a step-function forcing of wind stress, 
surface temperature, etc., produces results which imply that the eddy viscos- 
ity or diffusivity increases with time and tends to an infinite value for infinite 
time. However, this behaviour does not make the model invalid. We are here 
assuming infinite depth and neglecting the influence of stratification and ro- 
tation, and in such there is no natural limit to the size of turbulent 
eddies which may be produced. Such eddies may grow indefinitely, and the 
effective eddy viscosity or diffusivity has a term which increases linearly with 
time. If the eddy viscosity may be regarded as proportional to the product of 
a velocity scale and a length scale (the size of turbulent eddies), this will 
mean that such eddies should increase in size linearly with time, consistent 
with the lack of a natural limit to eddy size. For a constant wind speed U and 
surface bulk temperature T(0~), a steady state would have a constant current 
XU and a depth-independent temperature T(0~). However, this steady state 
will be reached very slowly in the assumed case of infinite water depth and 
zero Coriolis force. 



In addition, it is by no means assured that the vertical transport may be 
parameterised by means of an eddy viscosity or diffusivity. The model may be 
criticised for the fact that the velocity, temperature, or concentration gradient 
decreases to zero with increasing time: however, it is well known that it is not 
uncommon in turbulent flows for counter-gradient transport to occur, so that a 
finite flux of momentum/heat/mass under zero-gradient conditions should not 
be an unusual effect. It is common in the modelling of atmo spheric turb ulent 
diffusion for non-local transport algorithms to be used (e.g. IStulllll993h . and 



such algorithms almost invariably lead to counter-gradient transport in at 
least part of the domain, in agreement with results using large-eddy simulation 
algorithms ( Skvllinestad and Denbol 1995| ) . 



It is in fact possible to re-write Eq. 3 as an advection-diffusion equation, 1 



du du d ( du\ 
dt dz dz \ dz I 



employing the same value for eddy viscosity, ve = —v%z/(\U), as mentioned 
previously. The turbulent flux of momentum is thus expressed here as the sum 
of a downgradient and a no-gradient flux. The no-gradient flux is described 

1 H. Burchard, personal communication. 
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Table 1 



SkinDeEP deployment information for the measurements presented in this paper 



Date 


Stat. 


Station 


Times 


No. of 


Position 




no. 


name 


LST 


profiles 




1999-10-04 


4 


Punta Magdalena 


10:59-13:24 


120 


25°09.49'N, 112°59.52'W 


1999-10-10 


10 


T/S Irwin 


11:12-14:13 


161 


22°31.48'N, 109°35.43'W 


1999-10-14 


14 


Isla San Esteban 2 


06:15-07:44 


78 


28°36.15'N, 112°32.04'W 



by a downward vertical advection velocity w = —v%/(\U). Note that this 
advection velocity is in general much smaller than «*, so in our opinion it 
should not be regarded as unphysical. 



The eddy viscosity may still be described as the product of a velocity scale 
and a turbulent length scale —u^,z/(XU) = — A -1 (Cop a / Pw) 1 ^ 2 z. If we let 
the eddy viscosity profile be described in terms of a von Karman 'constant' 
K, i.e. v E = —ku*z, we have k = u*/(\U) = A -1 {CnPa/ Pw) 1 ^ 2 ■ This value of 
K is generally smaller than the 'classical' value of ~ 0.4: for the parameter 
values used here (see section 2.2), k increases with wind speed, being 0.050 at 
U = Om/s and 0.081 at U = 20m/s. 



3 Measurements 



A detailed description of the SkinDeEP instrument, which uses an FP07 ther- 
mi stor and a high-resolution plat inum wire (Pt) tem perature sensor, is given 
bv lWard and Minnettl l|2002h and lWard etafl (|2004bl ). The instrument is able 
to rise and sink autonomously, changing its density by inflating and deflating 
an external neoprene sleeve, and measurements are made during its ascending 
phase only. 



The instrument was deployed during the M OCE-5 cruise, in the waters around 
Baja California ([Ward and Minnettl 12004). Data were acquired from ten sta- 
tions: Table 1 summarises the deployment information for the three stations 
where the results are shown in this report. The measurements reported here 
are from the FP07 thermistor, as data from the Pt sensor were not available 
for this cruise. 



Other observations made include the following: 



• Sea-surface skin temperature by the M-AERI passive infrared radiometric 
interferometer, using the 500 -3000 cm -1 wavelength range, which has an ac - 
curacy of better than 0.05° dMinnett et al.ll200lUMinnett and Wardll2000h : 
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• Bulk water temperature measurements from the ship intake and from a 
floating sensor (an inverted hard plastic helmet filled with foam with a 
thermistor just below the waterline); 

• Air temperature, wind speed, and net heat flux. 

All three measurement series presented here are for periods with rather light 
winds, between 0.7 and 4.5 ms" 1 . The air temperature is less than the sea 
surface temperature (measured by SkinDeEP) on October 4, is greater than 
the surface temperature on October 14, and on October 10 the air and sea 
surface temperatures are quite similar. The net heat flux is positive (from air 
to sea) in all three cases. The net heat flux is primarily incident short-wave 
radiation, and is mo stly absorbed in the upper 0.1m of the water column 



fjSchlussel et al.lll990h . 



Shipboard logistics dictated that SkinDeEP was deployed during the afternoon 
to facilitate other operations, which is the timeframe when measurements 
were made on October 4 and 10. An exception to this routine was made on 
one day (October 14), when the profiler was deployed early in the morning. 
As measurements began shortly before sunrise, this was the only nighttime 
dataset, and was too limited to provide any conclusive evidence for daytime- 
nighttime differences. 



4 Model results 



Figures 3-5 show SkinDeEP measurements from three of the MOCE-5 sta- 
tions, together with background atmospheric parameters and net heat flux. 
Results of model simulations from equation 5 are shown at the bottom of each 
figure. The model is started with a uniform temperature at the time corre- 
sponding to the left-hand side of the plot, and the surface boundary condition 
T(0~,t) is set equal to the uppermost water temperature measured by the 
SkinDeEP profiler (at a depth of approximately 0.5 cm, in the bulk, below the 

1 /2 

surface skin layer). The value u* = (CnP n/Pv,) U in (5) is computed using 
the drag coefficient formula of[\¥y| JSJ: C D = A + BU, A = 0.8 x 10" 3 , and 



B = 0.065 x lO" 3 m^s. 

The ocean is assumed to be horizontally homogeneous, and neither the effect 
of stratification nor of the Coriolis force are taken into account. The current 
and waves were not calculated in this case, since no observations of current or 
waves were available. 

1999 October 4: 

In this situation the air temperature is lower than the sea temperature, but 
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O 22 -\ 
a> 21 




deg C 



22.4 
22.2 
22.0 
21.8 



21.0 



19 20 
1999 Oct 4 (UTC) 



Fig. 3. Observations and model results from Station Punta Magdalena, 1999 Oc- 
tober 4. 

the net positive heat flux should tend to increase the bulk temperature near 
the sea surface. The observed warming of the surface layer at around 12:30 is 
likely to be due to horizontal advection. Nevertheless, the model does predict 
an increase in depth of the warm surface layer (T > 22°C) to 2 m by the end of 
the observation period, which is of the same order as the depth increase which 
is actually observed. It may also be the case that the model parameterization, 
which predicts a 'von Karman coefficient' lower than the classical value of 
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19.5 



20 



20.5 



21 



1999 Oct10(UTC) 

Fig. 4. Observations and model results from Station T/S Irwin, 1999 October 
10, including details of the downwelling shortwave radiation (bottom frame). The 
modelled net heat flux is obtained from Eq. 6. The modelled bulk— skin temperature 
difference and CO2 flux are computed from the observed meteorological parameters 
from the ship (wind speed, air temperature, humidity, heat fluxes, precipitation) 
and uppermost water temperature measured by SkinDeEP, using the MOBLAM 
model of P. Schliissel and A. V. Soloviev. 
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Fig. 5. Observations and model results from Station Isla San Esteban 2, 1999 
October 14 

0.4, is underpredicting the amount of vertical mixing. In such further 
investigations are required to determine how the model parameters may be 
tuned. 

1999 October 10: 

Here the modelled increase in the warm surface layer depth during the ob- 
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servation period is quite similar to that actually observed. The drop in the 
observed surface temperature at 19:00 may be associated with the decrease 
in net heat flux in the previous 15 minutes, but the associated changes in the 
depth of the warm surface layer are not reproduced in the model. The greater 
wind speed in the period around 20:00, together with the positive air— sea 
temperature difference, tend to increase the temperature in the upper few 
centimetres of the ocean, and this is reproduced in both measurements and 
model results. Again, the agreement between model results and measurements 
is good, considering the fact that the model does not take stratification into 
account. 

The sudden cooling in the water surface at about 19:00 on October 10 coin- 
cided with the passage of a cloud, which is visible in the downwelling shortwave 
data (bottom frame of Figure 4). This had a dramatic and immediate effect 
on the surface warming, presumably because of the absence of wind mixing. 

1999 October 14: 

In this case the model appears to make a poor prediction of the temperature 
evolution of the surface layer of the ocean, when compared with the mea- 
surements. However, the temperature in the water column is rather uniform, 
varying by less than 0.15 K over the whole domain. The turbulent diffusion 
predicted by the model suppresses any developing non-uniformity in the tem- 
perature profile, an effect which would be mitigated if stratification were taken 
into account. 



Additional measurements and model results 

For the October 10 data set, a number of additional measurements and model 
results are presented. The downward heat flux Q m (t) through the surface 
computed from the present model is given by: 



with T(0 , t') being the uppermost water temperature measured by SkinDeEP. 

The bulk— skin temperature difference and gas flux are computed from the ob- 
served meteorological parameters from the ship (wind speed, air temperature, 
humidity, heat fluxes, precipitation) and uppermost (bulk) water temperature 
measured by SkinDeEP, using the Molecular Oceanic Boundary LAyer Model 




o 



f dT(0-,t') [ MQ) 2 
1 dt' [ \U{f) 



dt', 



(6) 
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(MOBLAM) 2 , by Peter Schliissel and A. V. Soloviev. The MOBLAM model 
is a surface renewal model, which takes account of solar and long wave radi- 
ation, turbulent and diffu sive fluxes of sensible and la t ent heat, and impact 
of rain on the cool skin ( Soloviev and Schliissell fl996L ISchfiissel et al. 1997; 



Craeve and Schliissel 1998). 



Heat flux: 



The heat flux computed by the present model, after a short period of adjust- 
ment, settles down to a fairly constant level which is roughly equal to the 
net heat flux deduced from observations. Given the simplicity of the model 
and the fact that no parameter adjustment has been made, this is remarkably 
good agreement. 



Bulk— skin temperature difference: 



The bulk— skin temperature difference calculated from the MOBLAM model 
gives values of between — 0.2K and — 0.15K, which are considerably smaller 
than the measured values using SkinDeEP and the M-AERI interferometer. 
This considerable discrepancy requires further investigation. 



Gas flux: 



The CO2 gas transfer velocity, computed by MOBLAM from the observed me- 
teorological parameters from the ship (wind speed, air temperature, humidity, 
heat fluxes, precipitation) and uppermost (bulk) water temperature measured 
by SkinDeEP, is fairly constant, except for enhanced values between 19:20 and 
20:00. The increasing values are associated by an increase in wind stress, which 
will increase the amount of turbulent diffusion, whereas the decreasing values 
are associated with an increase in the sea surface temperature and consequent 
decrease in CO2 solubility. 



Daytime-night-time differences: 



Shipboard logistics dictated that SkinDeEP was deployed during the after- 
noon. The one exception was on October 14, but the data set was too limited 
to provide any conclusive evidence for daytime-night-time differences. 



2 MOBLAM is available on the Internet from 

http://www.nova.edu/ocean/gasex/Moblam8.f90, the version used 

here being from 1999 February 10. The model is described on 
: / / www . nova . edu/ ocean/gasex/micro . html 
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5 Discussion and conclusion 



The time-dependent model presented in this paper is designed to provide an 
economical method of computing the evolution of the current (momentum), 
temperature, heat and mass flux in the near-surface layer of the ocean. Al- 
though it does not take account of stratification or of the Earth's rotation, 
nevertheless, given the simplicity of the model and the fact that no tuning 
was performed on the model parameters, it is shown to provide useful results 
for the time evolution of the temperature distribution in the oceanic surface 
layer over the measurement periods of the order of 2 hours. Some of the dis- 
crepancies we found may be explained as the result of horizontal advection, 
and by the model not taking account of the effect of stratification. It may also 
be the case that the amount of vertical mixing is underestimated, as indicated 
to a certain extent by the observations of Figure 3, and by the low von Karman 
'constant' of the model, discussed in section 2.4. It should be noted, however, 
that the turbulent flow and mixing in the near-surface boundary layer of the 
ocean need not be represented by the 'classical' von Karman 'constant' value 
of 0.4, as the physical conditions are by no means the same as those in a 
classical turbulent wall layer. A resolution of the problem of turbulent mixing 
intensity may be obtainable in the future by conducting an analysis of further, 
more detailed time-dependent coupled observations of velocity, temperature, 
and turbulence in the near-surface atmospheric and oceanic boundary layers. 

The computed net surface heat flux values are largely consistent with observed 
values. A suitable application for the model would be as an interface between 
complex time-independent models for interfacial flux of momentum, heat, and 
mass, and time-dependent three-dimensional numerical models of ocean circu- 
lation which employ turbulence closure schemes and account for stratification 
and the Coriolis force. 

The flux of CO2 and other gas species through the sea surface can be computed 
using MOBLAM or similar time- independent surface models. A straightforward 
extension of the present time- dependent model should be capable of extend- 
ing these flux predictions and computing the concentration and flux of CO2 
further into the water column. 

It should be possible to test the present model against time-dependent mod- 
els with on e (vert ical) spatial dimension, such as the warm layer model of 
Price et al. (|l986l ). and wit h models which pre dict bulk fluxes such as the 



TOGA-COARE algorithm (jFairall et al.lll996aU bh 3 . Such a study, including 



the use of turbulent closure models for stratified flows, and a more detailed 
evaluation of the effect of horizontal advection, would be a suitable topic for 



3 The algorithm is currently available on the Internet at URL 
http : //www. coaps . f su. edu/COARE/f lux_algor/f lux .html 
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a future paper. 
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